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USING AN IMPROVED HSDT-DEFORMATION THEORY IN 
ORDER TO BUILD UP A MATHEMATICAL MODEL FOR 
THE VIBRATIONS OF AN ORTHOTROPIC COMPOSITE 

BAR MAKING A SPATIAL MOTION 


Sabin RIZESCU', Dumitru EBANCA?, Dumitru BOLCU’, 
Marcela URSACHE*, Cora MIREA® 


Rezumat. Acest articol arata cum se poate construi un model matematic pentru 
vibratiile unei bare compozite avand simetrie elastica (ortotropica) considerand mai 
intai migcarea sa ca un solid rigid si, apoi, considerand o teorie a deformatiilor 
HSDT de ordinul al treilea in acord cu conditiile de compatibilitate Saint-Venant si, 
intr-un mod foarte original, a conditiilor Gay pentru miscarea barei ca un corp 
deformabil. In final este ardtat un algoritm bazat pe diferente divizate pentru a 
rezolva modelul. 


Abstract. This workpaper shows how to build up a mathematical model for the 
vibrations of a composite bar having elastic symmetry (orthotropic) considering first 
its motion as “rigid” body and, next, considering an HSDT-deformation theory of 
third order in full respect of Saint-Venant compatibility conditions and — in a very 
original way — of Gay conditions for the motion of the bar as a deformable body. 
Finally, a divided differences based algorithm designed to solve the mathematical 
model is shown. 
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1. Introduction 


There are a lot of attempts to model the elastic behavior and dynamics of 
composite materials [1], [2] and, especially the elastic behavior and dynamics of 
composite bars [3], [4]. 
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In [5] is shown how to build up a mathematical model for the vibrations of an 
orthotropic composite bar using the Hamilton’s variation principle and assuming 
some new and original hypothesis concerning the calculus of the elastic 
coefficients and the elastic displacements in full respect of the well-known Saint- 
Venant compatibility conditions of the Gay conditions [6]. We’ll show here how 
these results are to be applied in case of a composite bar making a particular 
motion. We’ll also show how the mathematical model obtained in this specific 
case could be solved using a divided differences based algorythm. 


A special focus is made on the motion of the bar as “rigid” body that turns out to 
be very important in order to establish the boundary conditions for the 
mathematical model. 


2. The motion of the bar like “rigid” body 


Let’s consider a kinematic linkage having two bars as elements (fig.1). The 
couplings from Oy and O3 are considered to be spherical. 


The linkage is rotating around the O,x axis with a constant angular velocity 


w=104,72 rad/sec. The reference system O,x/x$x° is accelerated on the 


> —>0 
vertical direction a, =—10,966 sin( 104,72t ) i, mis’, corresponding to a vertical 
>0 


motion described by the following equation: O,O, =0,001 sin( 104,72t )i; m. 


What really concerns is kinematics of the element No.2, considered as being a 
deformable one. We’ll show [7] that a motion that is analog to that when the 
element No.2 is considered as being one and the same with the element No.1 is 
fully compatible with dynamics of the linkage (considered for now as having rigid 
elements) despite the fact that couplings from O, and O, are spherical, so they 


are not introducing reaction moments. 

The motion of the element No.1 is considered as being imposed the way that: 
a=at; (1) 

where: ¢20 is time; 

@=const., @>0 is the angular velocity (its absolute value) of the element No.1 

= > 
Ca Sa): 
The angle 6 has - obviously - a constant value. 


We have established [7], [8] the motion like “rigid body” for the element No.2. In 
order to consider it as being deformable and to describe its vibrations we‘ll use the 
mathematical model established [8]. 
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This model has the following form: 


er bee] ele he ee eet 
(2) 


Vio aQee] él fysTe ol elbe} bee tes) 


Fig. 1 
Obviously we have: 
O02 
"@, = @,=0 
O0«2 00«2 
o,= 0,=0 


We'll focus on a specific case of a right composite bar having the next elastic 
characteristics: 
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- length: L,=0,45 m 

- dimensions of the its section: B= H =0,01 m 

- elastic coefficients: 

Eigg=2-7 410" Nim 
Eig = Ey = 3,18-10? N/m’? 
E4313 = E,33,=3,18-10° N/m? 

- the bar is considered to be an orthotropic one (even a transverse 
isotropic in order to make things more simple) and having elastic 
homogeneity: 

- the bar presents also mass homogeneity having the specific mass: 

p=1188 N/m’. 


We’ll consider the case of a velocity for the driving element of 1000 rot/min the 
way that: 


oA = 104,72 rad/ sec. 


Also we consider as known: 
- length of the driving element: L, =0,225 m 
- components of the acceleration of the origin O, with respect T, 
written in the T, reference system: 


22 2 2 
0 a, L, mr A Ay 


oa ie L > Fo,(2)2 =0; 
2 2 


o by o Leer a,L, 
ao,(2)3 = L iz 9 
2 2, 


- components of the specific mass gravitation force distributed on 


length: 
«| Damhe L 
oS 0 = 01188 [N/m] 


2 2 
3. Study of the vibrations of the composite bar 


exprimed in m/s” 


p, = -0,1188 


We are now ready to consider the bar as being deformable. We’ll consider the 
following initial conditions [8]: 


0 


Oy age ¢ e[07 |: {| 


0 
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0 0 
= 07; x, e|0,L,]; ‘a = 07; x, e|0,L,]. 
0 0 


t=0 


And the following boundary conditions: 


Oh, 


"|x, =O 


= (0) 


Cass) = 0; Chl Oe 6, (LSA 


; (V)t=0 
3x1 (We 


Having one single bar we’ll drop down the "j" index: 


Pf ea} feo la+ [Po +LeTtu}, 19 


(3) 
V Hols Yobebe lob} +b sto, {eh 
and, of course: x! =x, =x the way that we have: 
tu} = tu} (x1); {0} = {0} (x2) 
xe|0,L,|: te|o,r, | (4) 


The nature of the motion like “rigid body” gives {f} and {g} the following forms 
[9]: 


Uf} =Af} Gur); (gh = {const} (5) 
The matrix coefficients of the unknowns {u} and {a} are fully described in [8]. 


4. A divided differences algorithm designed to solve the mathematical model 


We are looking at the relation (3). We’ll divide [9] the interval [0, L2] into 
n subintervals having the pace on the h axis the way that a x, coordinate will be: 


x, =th; h=—; 0<i<n; (6) 
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We'll divide [9] the time interval of our study lo,r f |. t, being a final value, into 


m subintervals the way that at a certain moment of time f;: 


t 
Bie ‘hd ‘ 
t,=jt;T=—;O0Sj<m. (7) 
m 


The procedure (central differences) leads firstly to the following expression: 


P= He sard-2 bah Wi Help} fe} es P+ LP Lhe, H+ 
+P Heels tha te} eit} ® 


and: 
v= He, 218} 8, asst ile : 18. us} 0,, aS 
+ Teh fea ede Ee 
Wij) 218) +101) JA). (9) 
| u, =U (x,,t,); 0, aa bx.0)) fi =f (x,.2,),{e}= {const., 
Or, more: _ “Kh 


Seb Zell -Se kA PI. See) 
PT PT ae Te) SIP e.,,}+ 
Flee. nh ti} 
bb Lv ]|e Hak EV Th, {20'HV's 
+206) Ie..)-2 0 Ve.))-2Y Iu, }4 

+l fee: (10) 
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We'll note: 


Le} Hh 
Lp pL He) 


Sle} He. (11) 
and: 


Aye Hel 


ab hy} Lfy'|-[0)) 


2 VL WH} an 
With (12) and (11) the system (10) becomes: 


u an f= LAT" [B, le, tal’ [B, Jee, } s [4]’ [P*| Wier sf 


_f 


h? 
6 j}=[cT' by He, .}+lel'b. He,}- Stel] @.,}- 
~ssle' Wb) 8. }- lel bends ler Tt, }+[eP' te} 
1<i<n—1;j>0. (13) 
We have seen that: 
us = {Uj iUziU, } 
10} =10,50250,} (14) 


The boundary and the initial conditions leads us to the following conclusions: 
Wipf = (0},(V)i; (15) 
Wo, f= O}, (W) J (16) 
Wy S= 10}, (V) Js (17) 


= te ai fo tu, j-l 1, ~ {0} (18) 


So that: 
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Wir S= War fs (WD i 

0 }= {0}. (V)i 

I 

Fle) y= 0 
And more: 10. \ 10, \ (v) i 

1 

Os iF 0, j i =a = {0} 
And more 0), =9);5 (v) if 
in (46) having, obviously, ascendant differences. 

OAX,.t ;, =O 

“( (<i )-9, (x,t, )J= 0 

And, still more 0, (x, Ae ) oa (x, pet; ) 


Haslsyt,)-0: ly nt) =0 


or: 0, (eae )= 0; ke ..t,) 
in (24) and (23) utilizing descending differences. 


Let’s put in (13) j= i ei taking into account (18), (15), (20), (19) we have: 


}=[aP’[B, us} + TAT’ {fi} 
{9,}= cl'lb, |e, \ Ick fol: 1<i<n—I 
And because t=T: 
{u,}= Hens Ay, ITY th} 
A= [lle [p, | [c]’ {ge}; 1<i<n-1 


It becomes useful to make j =/ i into account (19) and (15): 


Win} = el : Ju - lA AJ '[P | UAT Be [P*| 4,1} 


= ch Oi} 


TP ]ie t+ Al t43 


YA 
mn 
Ly 
~zrlc 


Astra — j}+ [cl fe}s tsisn-1 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


nei h? tcl’) @..,,}- 


(27) 


Using an Improved HSDT-Deformation Theory in order to Build up a Mathematical Model 
for the Vibrations of an Orthotropic Composite Bar Making a Spatial Motion 119 


We'll note: bw, js 0: Gis Ca ay } : 
{w,,}= p: Oy n-1 p> 5,1 : } ; j=2,m-1 (28) 
Taking into account (28), (21) and (15): 


fu,.}=[AT"[B, Hee, }- SAPP? feo, }- 2 al P f+ 


he 2h 


2h 
112 } y [c}’ [D, ine - [c}’ IV" 0, J- “. [cy’ lv’ Ke, j- 
- lel] fai} [c}'{g}; 2<i<n—2 (29) 
The next values {u,,} and {6,,} continuing the algorithm as (27) imposes. 

With (17): 
en 2 } =[a]" IB. Jen 1 J- = [4y" [P*] {n> I }- = [a] [P’ Iov, J+ 

+s LAP PV, + FAD ihn) 
(6..2}=[cl' Bb. Ne..}-slct' Vn Folch] 6.2 

lel VT b2 leh te}: G0) 


And for 2 <i<n—2 we have at our disposal (27) putting j =/: 
On the other hand for i = 7 the system (13) becomes with (21) and (16): 
tH, ie } =[a]’ [B, trp }+ [a] [B, Ne, J- = A [P*| >, Js 
1 " I Y < 
~5, [4] ‘Pe. +1 ‘p+ |i6,,}+[4] Vir} 
0, itl j= [cy ID, ip J+ [ey [p, le, Js —|c}’ lV] 0, J- 
Lich le,,}-2 ich The, + IcP'ts an 
For 2<i<n-—2 the relations (13) remain available in the algorithm. 
For i=n-—Z, with (28) and (17) we have: 


bis es TAT BB, J fey rT Bs eo} SAP] fe} 


Sabin Rizescu, Dumitru Ebanca, 
120 Dumitru Bolcu, Marcela Ursache, Cora Mirea 


-- fall} SP 62 AT oe} 


2h 


(0, j+l } = [cy’ [D, (0... j-i he [c}’ [D, 0... j \. “ [cy’ aa bw, }_ 


Lich WIbe,s Zr hes Jel} 


for 2<i<n-—2 and using the system (13). 
This way we can determine the values for lus} and (0, \ 


The algorithm fits nicely to work it using programs like TurboPascal or C++. 


Conclusions 


We have shown how to build up and to solve a mathematical model for the 
vibrations of an orthotropic composite bar. More specific, the mathematical model 
will provide the evolutions in function of time for the displacement of the current 
point of the medium fiber “u” and for the rotation of the current right section “0”. 
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